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Abstract 

Entanglement of bosonic modes of material oscillators is studied in the context of two bi- 
linearly coupled, nonlinear oscillators. These oscillators are realizable in the vibrational-cum- 
bending motions of C-H bonds in dihalomethanes. The bilinear coupling gives rise to invariant 
subspaces in the Hilbert space of the two oscillators. The dynamics of entanglement and quadra- 
ture fluctuations is studied. Two classes of initial conditions, corresponding to the case where 
the total energy is concentrated in one of the modes and the other wherein the modes share 
the total energy, in the invariant subspaces are considered. The inadequacy of the known en- 
tanglement detection criteria is established and an inseparability criterion that is applicable to 
the states in the invariant subspaces is provided. Possibility of generating maximally entangled 
states is indicated. 



PACS: 03.67.Bg, 03.65.Ca 
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1 Introduction 



System with more than one degree of freedom has the potential to get entangled, a feature that 
is essentially quantal. No classical process has all the features exhibited by the correlations in the 
entangled quantum states. Such non-classical correlations are essential for quantum teleportation, 
quantum computation, etc. Hence it is important to identify physical systems wherein quantum 
entanglement is easily generated, controlled and measured. A host of proposals and demonstra- 
tions based on NMR, ion traps, SQUID, photon polarization, etc are known to generate entangled 
states. These systems can be used as gadgets to perform quantum computation [Il[2]. Recently, 
the entangling capabilities of molecules and the possibility of engineering their evolution to make 
quantum gates have been investigated [3l d]. Molecules have vibrational, rotational and electronic 
degrees of freedom. Often these degrees of freedom are coupled. It is precisely the coupling among 
the various degrees of freedom that is exploited to generate suitable molecular states. Study of such 
systems is important as the vibrational modes perform better as controllable qubits[4]. In addition 
to the aforementioned degrees of freedom, a nonlinear molecule like the dihalomethane has bending 
motion. Dihalomethanes are obtained by replacing two of the hydrogen atoms in CH4 with halogen 
atoms; CH2CI2 being an example. These molecules have nonplanar geometrical arrangement of 
atoms. A classical model to understand the vibrational spectra of these molecules takes the two 
C-H bonds to be two coupled quartic oscillators. The potential energy of the quartic oscillator is 

V{x) = ^Mw^rc^ - ax'^, (1) 

where x measures the deviation of C-H bond length from its equilibrium value. Here M is the 
"effective mass" of the C-H oscillator. The coefficient (a) of the quartic term is positive. The 
negative sign for the quartic term is to account for the dissociation of sufficiently excited molecule. 
Stretching of a C-H bond affects the charge distribution in the molecule. This, in turn, alters the 
potential of the other C-H oscillator, leading to an interaction between the two C-H bonds. Thus, 
the coupling between the bonds is the net effect of mechanical and electromagnetic interactions. A 
first order perturbation theory is performed to account for the interaction between the bonds. The 
small parameter e required for the perturbation theory is yJm/M^ where m is the electron mass. The 
mass of the hydrogen enters the analysis as the stretching of C-H bond stretching involves movement 
of hydrogen. A detailed discussion of these aspects and a derivation of Hamiltonian operator are 
given in [5|. The quantum Hamiltonian, obtained after making rotating wave approximation is, 



H 



[uj - '^){a)a + b^b) - ^{a^aa^a + b^bb^b) - e{a^b + aP) 



(2) 

Ho + Hn + H,. (3) 



This Hamiltonian describes the low energy vibrational states of the molecule correctly. The con- 
vention of setting ?i = 1 in the expressions is followed in this work. In calculations the actual value 
of the Plank's constant is used. The Hamiltonian H describes a bipartite system of two coupled. 



2 



bosonic oscillators. The operators a) and a are the creation and annihilation operators for the 
vibrational mode corresponding to one of the C-H oscillators. These operators satisfy the bosonic 



commutation relation 



a, d) 



I. The corresponding operators for the other C-H oscillator are ¥ 

and h. The bending motion is modelled by the term nonlinear in a) a and o^b. The term with coef- 
ficient w — 7/2 is Hq, the free Hamiltonian. The nonlinear part is denoted by Hn and its coefficient 
is defined by 7 = 3a?i^/4M^a;^. The coupling term is Hi which has e = ^m/M as the coefficient. 
In the context of nonlinear optics this type of bilinear coupling has models nonlinear couplers for 
electromagnetic fields [H [7j. The Hamitlonian has been used in the study of local modes (which 
are the states of molecules where energy is concentrated in one of the C-H bonds in the molecule) 
and quantum lattice solitons[5]. In the present work, the Hamiltonian H is analyzed to study en- 
tanglement in the coupled bosonic system. Analysis of bosonic systems in the context of quantum 
information theory is becoming important[8], especially after the demonstration of teleportation of 
coherent states of electromagnetic fields[9]. Theoretical studies have indicated the possibility of gen- 
erating gates with bosonic modes and the robustness of these gates against thermal dissipation [TO] . 



The Hamiltonian in Eql2] is similar to that used in the analysis of propagation of two-mode 
electromagnetic fields in a medium with \^ nonlinearity [1 1 1 [T2] . The nonlinear medium couples the 
two modes of the electromagnetic field propagating in it [13]. An important difference between the 
Hamiltonian for a nonlinear coupler and H is that the coefficient of the nonlinear term is positive in 
the former case while it is negative in the later. The Hamiltonian has been investigated in the study 
of systems for generating maximally entangled states [14^ 115] , entanglement dynamics [16] , control 
of switching modes [T7] , generation of bright entangled continuous variable states [18] , wavepacket 
dynamics [19], nonlinear quantum scissors to generate finite dimensional states in systems with 
infinite dimensional Hilbert space [20jand many more interesting physical applications. 

In the case of dihalomethanes, the Hamilotnian describes mechanical oscillators, that is, oscil- 
lating masses. There are no external agents, like the Kerr medium in the case of nonlinear couplers, 
responsible for the coupling between the two oscillators. The coupling is internal and it is not pos- 
sible to switch off the coupling. The presence of coupling alters the dynamics of the two modes 
corresponding to the two C-H oscillators significantly. In the present work, the physical properties 
such as quadrature fluctuations and entanglement are studied. The studies pertain to the situation 
when dissipation and external fields are absent. Effects of the nonlinearity and the coupling on the 
dynamics are explored in this work. The results presented here correspond to typical values of 7, lv 
and e for a class of dihalomethanes. The relevant values are tabulated in Table. [I|5j. The values 
of the parameters are nearly equal for the three species. For the purpose of presenting the results, 
the representative values used for the parameters 7, e and u are 125, 30 and 3050 respectively, all 
expressed in units of cm~^. The organization of the paper is as follows. In Section [2] the spe- 
cial features of the Hamiltonian are discussed. In Section [3l dynamics of entanglement is studied. 
Known entanglement detection criteria are analyzed in the present context. The inadequacy of 
these criteria is established and a new criterion suitable for states of the present system is given. 
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Table 1: Values of 7, e and u expressed in cm units 



Molecule 


7 


e 




CCI2H2 


127.44 


29.54 


3020.1 


CBr2H2 


125.45 


32.80 


3026.8 


CI2H2 


124.25 


33.69 


3068.7 



Results on the fluctuations of single-mode and two-mode quadratures are presented in Section [H 
In Section [Sj it is shown that maximally entangled states are generated at specific instants during 
evolution from a separable state. Results are summarized in Section [6l 

2 Special features of H 

The presence of the bilinear coupling term He endows the Hamiltonian H with interesting features. 
A suitable basis to expand an arbitrary state of the oscillators is the set {|n, m)}, n, m = 0, 1, 2, • • •. 
The quantum numbers n and m label the states of the oscillators corresponding to the two modes. In 
the absence of the coupling, for n + m = N, then there are + 1 eigenstates of the form \N — m, m) , 
with m ranging from to N. The corresponding eigenvalues are {lo — ~ 2 [(^ ~ '^)^ ~^ ■ 
For odd A, there are (A + 1)/2 distinct eigenvalues and the states are doubly degenerate. When N 
is even, there are [N + 2)/2 distinct eigenvalues; there is one nondegenerate state \N/2,N/2) and 
the rest are doubly degenerate. When e 7^ 0, the degeneracy is lost and the states |n,m), except 
when n = m = 0, cease to be eigenstates of H. The free Hamiltonian Hq, which is the total num- 
ber operator, commutes with the total Hamiltonian H and the coupling term Hi. Consequently, 
the Hilbert space of the coupled system splits into disjoint, irreducible, invariant subspaces. Each 
invariant subspace is characterized by total the quantum number N = n + m and the symbol Sn 
is used to represent the corresponding invariant subspace. The subspace Sn is the span of the 
vectors {\N — m, m)}^^Q and its dimension is + 1. For instance, when n + m = 1, there are 
two possible states, namely, |0, 1) and |1,0). The span of these two canonical basis states forms 
the relevant invariant subspace ^i. If coupling is absent the irreducible, invariant subspaces are 
of dimension one. The nonlinear and the coupling terms in the Hamiltonian do not commute. 
Consequently, the factorization of the time-evolution operator is difficult. Nevertheless, including a 
nonlinear coupling term, namely, a)a¥b, allows to write the Hamiltonian as a sum of three mutually 
commuting terms [16]. In the present work the discussions are limited to the bilinear coupling model. 

Another remarkable feature of the invariant spaces is that the states of the form |n, m) are 
the only product states in the respective invariant subspaces. Every other state in a given invari- 
ant subspaces is entangled. The most general state in the invariant subspace Sn is of the form 
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J2n=o (^nln, N — n). If this is to be a product state, then it should be expressible as J2n=o Cn\n, N — 
= J2k=o 9k\k)a X Z]m=o /ml"^)f" suffixes in the states indicating the corresponding modes. 
Consequently, the coefficients satisfy c„ = QufN-n- Let the state of the bipartite system is such 
that at least two coefficients, say, Cp and Cq are nonzero. The relation among the coefficients im- 
plies that gp, Qg, f]\f-p and /n-q are nonzero. As a result, the product state is gpf]\f-p\p, N — p) + 
gqfN-g\q,N-q)+gpfN-g\p,N-q)+ggfN-p\q,N-p). However, the states \p,N-q) and \q,N-p) 
do not belong to Sjy, their total quantum numbers are N — q + p and N — p + q respectively and 
not A^. Hence, it is not possible to express the state with two nonzero c„ as a product of states 
in Sn- This argument can be extended to states wherein more number of c„ are nonzero to show 
that such states are entangled. If there is only one non-vanishing c„, then the gn and fN-n are the 
only nonzero coefficients. The corresponding state \n,N — n) is a product state. 

Explicit construction of the eigenfunctions in terms of the number states and eigenvalues of the 
Hamiltonian is possible by the method of number states [5]. However, the resulting expressions for 
the eigenstates are complicated. On treating the coupling as a perturbation, simple expressions 
for approximate eigenstates and eigenvalues are obtained by first order perturbation theory. Since 
{N — m, m\a^b + ab^N — m, m) is zero, the eigenvalues do not change to first order in perturbation. 
Apart from an overall nonramlization factor, the perturbed eigenstates, denoted with a suffix p, 
are 

I A — m, m)p oc \N — m, m) + /i| A — m -|- 1, m — 1) -|- /2I A — m — 1, m + 1) (4) 

where 



_ e a/(A - m + l)m 
~ 7 l + A-2m ' 



e v/(A-m)(m+I) 
- 7 l-A + 2m • 

If |A — 2m\ = 1, either fi or /2 becomes infinity and the perturbation theory is not applicable. 
Though the eigenstates are doubly degenerate, the results of nondegenerate perturbation theory 
are applicable if | A — 2m| 7^ 1. This holds for all for states except those states of the form |m, m+ 1) 
and \m + 1, m). The transition matrix element {m + 1, m\ab^ + a^b\m, m + 1) is nonzero and hence 
the nondegenerate perturbation expansion is invalid. The matrix element of -|- a^b among the 
other states vanish and the EqnH] holds for such states. Further it is required that 7 >> e for the 
perturbation expansion to be valid and this condition holds in the case of dihalomethanes. Two 
important special cases, |A, 0) and its degenerate counterpart |0,A), are analyzed. Physically, 
these states correspond to the situation in which the total energy is concentrated in one of the 
modes and are referred as "local modes". From the expression for perturbed states in Eqn. [U 



\N,0)p oc |Ar,0) + ^^^|A-l,I), (7) 

\0,N)p oc |0,A) + --^|1,A-1), (8) 
7 A — 1 
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where > 1 to ensure that the perturbation results are vahd. For large N, the second term of 
the perturbed state is negligible. Consequently, the unperturbed states |A^, 0) and |0, A^) qualify as 
approximate eigenstates of the perturbed Hamiltonian H. The validity of the result is ascertained by 
studying the overlap between the initial state and the corresponding time-evolved state. Fidelity 
is defined as the absolute value of overlap of the initial state with the time-evolved state. This 
definition means that the fidelity is given by the absolute value of the expectation value of the 
evolution operator exp{—itH) evaluated in the initial state. If a state is an eigenstate of H then its 
fidelity is unity at all times. The fidelity for the states |2, 0) and |1, 1) are given in given in Fig. la. 
The results are based on nonperturbative evolution of the initial state so that the applicability of the 
perturbation theory is justified. The explicit construction of the Hamiltonian matrix required for 
the numerical computation is specified in the next section. The fidelity of the state |2,0) becomes 
very small, an indication that the state is not an eigenstate. In the same figure the fidelity of the 
state |1, 1) is given. Though the fidelity does not become zero, its deviation from unity is significant 
as the state is not an eigenstate. In Fig. lb fidelity of the states |0,4) and |2,2) are given. It is 
seen that the fidelity of the former remains close to unity for all times. Thus the state |4, 0) is an 
approximate eigenstate. This is expected based on the perturbation theory result that for large 
the local modes are approximate eigenstates. The fidelity of the state |2, 2) is appreciably different 
from unity since the state is not a good approximation to an eigenstate. Though the state |1, 3) has 
significant amount of energy concentrated in one of the modes, the total quantum number (equal 
to 4 in this case) is not sufficient to make it an approximate eigenstate. 

The role of nonlinear term in the Hamiltonian deserves to be stressed. For the local modes 
to be approximate eigenstates the nonlinearity parameter 7 should be large in comparison to the 
coupling strength e. In addition, the eigenstates of H are all entangled, except the ground state, 
if nonlinearity is absent or weak. Nonlinearity makes it possible for the system to have approxi- 
mate eigenstates which are separable, despite the presence of coupling. If 7 vanishes, first order 
perturbation theory fails and the states of the form \n,m) do not approximate the eigenstates. 

3 Dynamics of Entanglement 

In the absence of coupling, the Hamiltonian is a function of the two number operators, namely, 
a^a and b'^b. The states — m,m) are eigenstates of the coupling-free Hamiltonian. From the 
expressions for the perturbed eigenstates given in Eq. SI it is clear that if the nonlinearity dominates 
over the coupling, 7 >> e, the separable states |iV — m,m) become the approximate eigenstates of 
the total Hamiltonian. The other extreme case corresponds to the absence of the nonlinear term 
while the coupling is retained. Perturbation theory results are not valid in this case. However, in this 
limit the Hamiltonian H is the sum two operators, namely, a)a-\-o^b and a)b + ab^ , which commute. 
By the Baker-Campbell-Hausdorff formulap^. the unitary time evolution operator exp(—itH) 
factorises to 

U = exj){-itH) = exp{iet{a^b + ab^)) exj){-iujt{a^ a + b^b)). (9) 
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10 20 30 40 50 60 70 80 90 100 10 20 30 40 50 60 70 80 90 100 

Time (ps) Time (ps) 

(a) N=2 (b) N=4 

Figure 1: Fidelity {exp{—itH)) of states corresponding to total quantum numbers N = 2 and 
= 4. The expectation value is cvahiatcd in the states for which fidelity is calculated. Fidelity of 
|0, 2) (dashed) and |1, 1) (continuous) are shown in (a). Fidelity of |0, 4) (dashed), |2, 2) (continuous) 
and |1, 3) (dotted) are shown in (b). Values used for the parameters are 7 = 125 cm~^, e = 30 cm"-*^ 
and Lo = 3050 cm^^. 
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The time-evolved states are obtained by the action of U on the initial states. 

For dihalomethanes, the parameters 7 and e satisfy 7 ~ 4e (see Tabled]). Hence, the advantages 
of the limiting cases, namely, 7 = or e = 0, in solving for the dynamics are not available. However, 
the availability of invariant subspaces simplifies the study of dynamics if the initial state is in one 

the evolved state has no overlap 



of the invariant subspaces. If the initial state is ^2^=0 Crl-^ 
with states in other invariant spaces. Hence, the dynamics of such states is dictated by a truncated 
Hamiltonian of dimension + where N is the total quantum number characterizing the invariant 
space containing the initial state. The Hamiltonian, expressed in the canonical basis \n,m), is the 
sum of a diagonal matrix and a non-diagonal matrix. The diagonal matrix is the sum of the free 
Hamiltonian Hq and the nonlinear term Hn- The corresponding matrix elements are 



■r, r 



(Ho + H, 



n )j,k 



2Nj + 2f) 



Oj,k 



j,k = 0,1,2,- --N. 



(10) 



The interaction He is nondiagonal whose elements are 



(FeWi = v/(i + i)(^-i) j 



1,2,- • -N-l, 



(11) 



and other elements vanish. This truncated Hamiltonian matrix is sufficient to study the evolution 
of states belonging to the particular invariant subspace S'at. This finite size matrix makes it easier 
to numerically track the evolution of the system. The results given subsequently are based on 
numerical computations. For the calculations presented here the largest matrix is of order five, 
used to study states in the invariant subspace 5*4. For an initial condition involving states from 
two or more invariant subspaces, the Hamiltonian matrix has to be enlarged. It is to be noted 
that the evolution dictated by H is solvable analytically [25]. If a{t) and b(t) represent time-evolved 
operators, they are related to the initial operators a and b through 



^Hi)t] + H.c. 



exp[— i(e 

exp[—i{e — — H.c. exp[— i(e 



exp[-z(e - ^Hi)t] 



H.c. 



fHi)t] + H.c. 




(12) 



Here Hi = -|- a6^. The term exp{i^tHi/2) does not factorize further. As a consequence, the 
expressions for expectation values become lengthy and cumbersome. This justifies the use of nu- 
merical calculations which are very accurate. 



As noted earlier, states of the form \N — m, m) are the only product states in the invariant subspace 
Sn- Due to the coupling, the two modes may get entangled during evolution even though the initial 
state is a product state. Quantification of entanglement between the oscillators is easily done using 
density operators. If the initial state of the system is in one of the invariant subspaces, say, Sm, 
then the density operator for the bipartite system is 

N N 

p=^^Cncl^\N -n,n){N -m,m\. (13) 

m=0 71=0 
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The coefficients are time-dependent and satisfy suitable initial conditions. The reduced density 
operators for the two oscillators are 



N 



Pa = Y. \cn?\N -n){N -n 



(14) 



n=0 



and 



TV 




(15) 



n=0 



The suffixes a and h indicate the two modes respectively. The reduced density operators are 
diagonal in their respective number state basis. Further, the diagonal elements of the reduced 
density operator are the same except for a reversal of ordering: the coefficient |cnp appears as 
the probability for the state \n) for the a-mode while it appears as that of the state \N — n) for 
the other mode. This restriction is only because the state of the bipartite system belongs to the 
subspace Sn- 

A measure of entanglement in bipartite, pure states is the linear entropy L defined as follows [21]: 



The numerical value is independent of whether pa or is used in the expression. Using the reduced 
density matrix expression given in Eq. [T^ the linear entropy is 1 — Y^n=o |cn|^- Another measure 
of entanglement [2T] is von Neuman entropy S defined as 



Its value is zero for separable states and reaches the maximum value of log2 M for a system whose 
Hilbert space is of dimension M. The explicit expression for S is 2^)^=0 I log2 |c„,|. Nonzero 
values of entropy {S or L) imply that the system is entangled. In Fig. 2 the dynamics of S is 
shown as the system evolves from initial conditions of the form |n,m) which are separable states 
and hence the entropy is zero. Instead of plotting 5, the ratio of S and log2(A^ + 1) is plotted. The 
ratio varies from zero to a maximum of unity for any state. This makes it easier to compare the 
evolution of S for initial states from different invariant subspaces. The evolution of entropy when 
the total quantum number is one, two and three are shown in Figs. 2a -2c respectively. The en- 
tropy increases substantially approaching the maximum attainable entanglement and subsequently 
oscillates. However, the state |0, 4), whose entropy evolution is shown in Fig. 2d, does not evolve 
to states of high entropy. The entropy of the local mode |0, 4) is zero as it is a product state. Being 
an approximate eigenstate, the evolved states do not differ significantly from the initial state. As a 
consequence, the entropy remains low during evolution. In practical terms, this implies that local 
modes are not good candidates for generating states of high entanglement using bilinear coupling. 
The situation changes drastically if other initial conditions are chosen. When the system evolves 



L = l- Tr [paPa] ■ 



(16) 



S = Tr [pa log Pa] . 



(17) 
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from initial conditions other than the local modes, the entropy S attains values closer to the maxi- 
mum allowed value of log2 N. Evolution of entanglement when the initial states are not local modes 
is shown in Fig. 3. In particular, comparison of Fig. 2d for the state|0,4) and Fig. 3(d) for the 
state 1 1,3) of same total quantum number, shows the clear distinction between the evolution of 
a local and a nonlocal mode for = 4. The evolution takes the system to maximally entangled 
states whenever the reduced density operator of the form J2n=o during evolution. 

It is important to note that if the nonlinear term is absent, the bilinear coupling cannot lead to 
entangled states from separable, initial states unless the initial states are nonclassical[22t . Here 
"nonclassical" implies that the Glauber-Sudarshan P- function for the state in the diagonal coherent 
state representation does not qualify as a probability density [Ml [26]. The states in Sn for N > 
are non-classical and hence the evolved states are entangled. The presence of the nonlinearity 
along with the bilinear coupling can entangle states which are classical. In essence, bilinear cross 
coupling, without the presence of nonlinear term in the Hamiltonian, cannot generate entanglement 
if the initial states are classical. 



0.8 

g:0.6 

o 

iS 0.4 
0.2 


20 40 60 80 100 20 40 60 80 100 

Time (ps) Time (ps) 





Figure 2: Temporal evolution of von Neuman entropy S and linear entropy L. The quantities plotted 
are '^og'^(JV+i)^ (continuous curve) and L (dashed curve) for various initial conditions. Figures (a)- 
(d) show the evolution as the bipartite system evolves from the initial conditions |0, 1), |0, 2), |0, 3) 
and |0, 4) respectively. Values of 7,e and uj are as indicated in Fig. 1. 



The states under consideration are non-Gaussian, pure states. While von Neumann entropy 
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Figure 3: Temporal evolution of von Neumann entropy S and linear entropy L. The quanti- 
ties plotted are '^og^(jv'^i)^ (continuous curve) and L (dashed curve) for various initial conditions. 
Figures (a)-(d) show the evolution as the bipartite system evolves from the initial conditions 

|1, 0), |1, 1), |1, 2) and |1, 3) respectively. For the state |1, 3). the ratio attains values close to 0.8, in- 
dicating the entanglement is nearly 3 ebits. Higher energy states can generate more entanglement. 
Values of 7, e and u are as indicated in Fig. 1. 
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provides the unique measure of entanglement for pure states, it is required to have a criterion 
expressed in terms of measurable quantities like the moments of the creation-, annihilation- and 
number- operators. Many criteria have been proposed to detect entanglement, but none of them 
are universally applicable. All these criteria are sufficient conditions for entanglement. For in- 
stance, the sufficient and necessary conditions for entanglement in bipartite Gaussian states are 
known [271 128j . However, for non-Gaussian states these criteria may not work. In what follows, 
several of the known criteria are applied to the states in the invariant subspaces. The explicit 
calculations are simplified by using the result that the expectation value of a^fe^atPftt? states 
in any invariant subspace is zero if m -\- n ^ p + q, where m, n,p and q are non-negative integers. 



1. Duan et al criterion [27]: This criterion is expressed in terms of uncertainties in the two-mode 
operators: 



u 



1 

1 

iV2 



|m|(a -|- dt) -|- 



m (a — a 



m 

t)-l(6-6t) 
m 



(18) 
(19) 



Separable states satisfy ((An)^) + ((Aii)^) > (1 -|- m^)/m^. For the states in Sn, the relation 
becomes ((An)^) -|- ((Aw)^) = (|m|^ + l)/m^ + \m\^Na + Nh/m?, where A''^ and refer to the 
expectation values {a) a) and {b'^b) respectively. Since the number operators are non-negative op- 
erators, the sum of the uncertainties indeed satisfies the condition of separability although the 
states are entangled as evident from the earlier discussion on the von Neumann entropy. Hence 
this criterion is not useful in the present context. 



2. Mancini et al criterion : This criterion employs the product of uncertainties in the operators 
considered in the criterion due to Duan et al, choosing m = 1 ^29j . The inequality to be satisfied to 
identify entangled states is that 

((A^)2)((At;)2) < 1. (20) 

When applied to the states in Sn, ((Au)^) = 1 + (d^d + wb) + (d^^ -|- ao^) which can be written as 
1 -|- {{a) + o^){a + b)). The operator (d^ -|- b^){a -\- b) is of the form A^A, which is positive. Hence, 
((Au)2) > 1. Similarly, {{Av)'^) = 1 + ((dt - U){a - b)) > 1. Thus, the product of the uncer- 
tainties is larger than unity making the criterion not useful in identifying the entangled states in Sn- 



3. Shchukin-Vogel's criteria: 

Many of the known criteria are expressible in terms of the determinants introduced in [30 
All separable states satisfy 

1 (dt) (6t) 
= (d) (dtd) (dtftt) > 0. 
(6) (dS) (6t6) 



(21) 



12 



This inequality is a stronger version of Duan et al criterion. If the determinant is negative, the state 
is entangled. For the states in an invariant subspace Sn, the determinant is explicitly evaluated to 
give {a^a){b^b) which is positive. Thus, the criterion does not detect entanglement in these states. 
Another criterion that works for entangled coherent states \a)\P) — \ — a) \ — P) is that 



1 (6) (a6t) 
(at6) (a^b^b) {a)ab%) 



<0, 



(22) 



implies entanglement. For the states under consideration in the present work, the determinant is 
{a^awb){wb) which is positive and hence inconclusive about the entanglement in the states. 

4. SU(2) criterion: The operators corresponding to the two modes can be combined to obey 
SU(2) and SU(1,1) algebra. The corresponding uncertainty relations provide criteria to detect en- 
tanglement. The expression for SU(2) uncertainty product satisfied by separable states [HI [32] 
when applied to the states in the invariant subspaces is 



{a^aV^b) [4(0+06+6) + 2(d+a) + 2(6+6)1 + 4(a+a)(6+6) > 0. 



(23) 



However, to detect entanglement this expression should attain negative values. Since the expression 
is positive, it does not detect entanglement in the states. 

5. It is known that SU(1,1) [311 [32] uncertainty relation can detect entanglement in a class of 
two-mode non-Gaussian states involving the ground and the first excited states of harmonic oscil- 
lator. For the states in any of the invariant subspaces, the criterion states that if 



/ll = [l + 2Nab + Na + Nk-2NaNk] 

-|iVa + 7Vb)|2 <0, 



3f?(d26+2)2 - 5R(a+6)2 



(24) 



the corresponding state is entangled. Here Nah is the expectation value {a)aWb) and 3f? stands for 
real part of the expression that follows. 

Simiarly, Simon's criterion [28j identifies entangled states in the invariant subspaces if 



fs 



^{ab^y 



+ \{abX - \{ab^)\^ [Na + iVfe + 2NaNt] 



1 



lb 



The criterion due to Hillery and Zubairy [33] states that if 

fHZ = NaNh-\{a^b)\^ <0, 



(25) 



(26) 
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then the state is entangled. For the states generated when the system evolves from the initial 
condition |2,2), these three criteria are evaluated numerically. The results are shown in FigUl It 
is seen that the entanglement is not identified by the criteria fnz and fs ■ The SU(1,1) criterion 
detects entanglement in at least some of the states occurring during evolution. This is the reason 
for choosing the |2,2) as the initial state. If any other product state in 5*4 is chosen as the initial 
condition, all the three criteria fail. 
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Figure 4: Temporal evolution of criteria for entanglement as the system evolves from the initial 
state 1 2, 2). The criteria shown are /i^s (dashed) , fs (dash-dot) and /n (thin continuous). For 
comparison, the von Neumann entropy is shown (thick continuous). The curves correspond to ratio 
of the criterion to its absolute maximum value during the evolution. This quantity has been chosen 
so that the various criteria are limited to a maximum value of one. This ratio is sufficient as it is 
the occurrence of negative value that is of importance and not actual values of the criteria. The 
value of /ii falls below zero (the horizontal line) between 800 ps to 1000 ps indicating entanglement 
during those instants. At other instants, the criterion does not detect entanglement. Values of 7, e 
and uj are as indicated in Fig. 1. 

The failure of the known criteria to detect entanglement provides motivation to look for a new 
criterion. The states in an invariant subspace are spanned by states of fixed total quantum number 
N = n + m, leading to strong correlation between the number of quanta of the first mode and 
that of the second mode. It is, therefore, natural to expect that a correlation function of the 
number operators of the two modes may detect entanglement. Such a criterion for entanglement 
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is derived as follows. For all product states the equality {a)ab^h) = {a)a){b^h) holds. Consider the 
most general state '^Cn\n,N — n), where the summing index n runs from zero to N. For these 
states, {a^ab^b) = N{a^a) - {a)aa)a) and {a^a){b^b) = N{a''a) - {a^a)'^. But, {a^aa^a) > {0)0)'^ for 
any state. Comparing the expressions for {a^a¥b) and (0^0) (6^6) yields the inequality {a^ab^b) < 
{a^a){¥b). The equality sign holds if (a^aa^a) = {0)0)"^, which is true if all but one of the c^s is zero. 
In that case the state is a product state and there is no entanglement. To bring out these features, 
the temporal variation of the number correlation function D = {a^aWb) — {a^a){b^b) is shown in 
Figini Negative value of D implies that the state is entangled. For comparison, the von Neumann 
entropy is shown. It is seen that whenever the entropy is positive, which indicates entanglement, 
the quantity D is negative. Thus, D qualifies as a suitable criterion to detect entanglement in the 
class of non-Gaussian, pure states considered. 



1 r 




Time (ps) 

Figure 5: Temporal evolution of criteria for entanglement as the system evolves from the initial 
state 1 2, 2). The quantities plotted are von Neumann entropy (conitnuous) and D (dashed). Values 
of 7, e and lo are as indicated in Fig. 1. 



4 Quadrature variances 

The coupled C-H oscillators modelled by H are oscillating masses. The quadratures of these 
oscillators are naturally identified with their positions and momenta. For the individual oscillators 
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described by H, the position and momentum quadratures are given by 

Qa 



Pa 



a + a) 
a — 

ITT' 



(27) 
(28) 



respectively and they satisfy [Qaj -fa] = where fi is set equal to unity. The corresponding variances 
satisfy 2A(5aAPa > 1. If the variance of one of the quadratures is less than 1/2, the corresponding 
quadrature is said to exhibit squeezing. Evaluation of the variances in the present case is simplified 
since the expectation values of operators involving creation and annihilation operators with unequal 
exponents, vanish in the states within an invariant subspace. For instance, the expectation value 
of ac^'^ is zero as creation and annihilation operators occur with different exponents; similarly the 
expectation values of ah^a^h\ etc all vanish for the states in Sj^. Hence, for the states belonging 
to Sat , 



and 



(29) 



(30) 



The suffix a labels the mode. Similar expressions can be written for the other mode and the results 
on variances are identical to those of the a— mode. If A'^ > the quadratures of the individual oscil- 
lators do not exhibit any squeezing for the states in the invariant subspaces. Further, the numerical 
values of the variances in the two quadratures are equal and increases with the total quantum num- 
ber N. The state |0,0) does not exhibit squeezing and it corresponds to minimum uncertainty 
state. The notion of quadratures associated with the individual modes has been generalized to 
multi-mode cases. For the bipartite system, the two-mode quadratures are defined as[ 



d2 



+ b + a^ + U 

3 : 
22 

+ 6 - gt - 6t 

3 • 

i22 



(31) 

(32) 
(33) 



Based on these expressions for the two-mode quadratures, the respective variances of the quadra- 
tures in any state belonging to Sn are given by 



(Ad2)' 



1 + N+ (afot) + {a)b) 



(34) 
(35) 
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The states belong to the invariant subspace Sn possess equal two-mode quadrature variances. The 
two-mode quadratures satisfy the commutation relation [1^1,^2] = i/2 and the corresponding un- 
certainty relation is AdiAd2 > 1/4. If any of the quadrature has uncertainty lower than ^, the 
quadrature is said to be squeezed. The fact that the two-mode variances are equal implies that one 
of the variances cannot be reduced to values less than 1/2, the condition for squeezing, without 
violating the uncertainty relation. In short, squeezed fluctuations in the single-mode or two-mode 
quadratures are not possible for the states in any of the invariant subspaces. Though there is no 
squeezing, the states in the invariant subspaces are very much nonclassical in the sense that their 
respective P-distributions are highly singular. 



5 Generation of maximally entangled states 

Bell states enjoy a special status in quantum information theory as they are maximally entangled 
bipartite states[T]. The canonical Bell states are the linear superpositions |0, 0) it |1, 1 and |0, 1) it 
|1, 0). Bell-like states are similar combinations allowing for a relative phase between the superposed 
states. Though the C-H oscillators are not two level systems, the presence of invariant subspaces 
allows for the creation of Bell-like states for the two coupled oscillators. This special feature is 
available in the invariant subspace 5i spanned by |0, 1) and |1,0). The Hamiltonian H2 required 
to describe the dynamics of the states in the subspace ^2 is 



H2 



w — 7 — e 
— e w — 7 



The eigenvalues of the H2 are a; — 7 it e and the corresponding eigenstates are the Bell states 
|±) = [|0,1)±|1,0)]/V2. 

If the initial state is |0, 1) = [|+) + |-)]/\/2 then 

exp{-itH2)\0,l) oc (1 + exp(iet))|0,l) + (1 - exp(iet))|l, 0), (36) 

apart from an overall multiplicative normalization factor. If the time t is chosen to be et = tt/2, 
then the time-evolved state is the Bell-like state |0, 1) 0), apart from an overall phase factor. If 
et = 37r/2, the initial state evolves to become another Bell-like state |0, l) + i|l,0), but for an overall 
phase factor. The overlap between the time-evolved state of the system and the Bell- like states are 
given in Fig. [6l At specific instants, the initial state evolves to have unit overlap with |0, 1) +i|l, 0); 
there is an overall phase factor to the actual state achieved and the overlap is insensitive to such 
factors. Since the subspace 5*1 is invariant under the unitary evolution, the other two Bell-like 
states, |0,0) it |1,1) are not attainable with any initial condition contained in 5*1. These states 
belong to the direct sum of Sq and 5*2. Without external control fields it is not feasible to generate 
all the Bell-like states using the inherent coupling in the system. 
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Figure 6: Overlap between the state exp(—itH2)\0, 1) and the Beh-hke states |0, 1) 0) (contin- 
uous) and |0, 1) — i|l,0) (dashed). The instants when the overlap becomes unity, the initial state 
has evolved into the corresponding Bell-like state. Values of 7, e and lo are the same as for Fig. 1 



18 



6 Summary 



The Hamiltonian for the coupled, nonhnear C-H oscihators of dihalomethanes is similar to that of 
the nonlinear couplers with negative x^. Thus the C-H oscillators are the microscopic nonlinear 
couplers for bosonic fields of vibrational-cum-bending motion. Like the Kerr nonlinear couplers, 
the dihalomethane can serve as a system wherein entanglement is easily generated with the advan- 
tage that there is no need to have external fields for the purpose of entangling the modes. The 
coupling due to the bending motion of the molecule creates entanglement between the modes and 
allows for exchange of energy. The structure of the Hamiltonian with commuting terms allows 
to split the Hilbert space into irreducible, invariant subspaces. Each of the invariant subspace is 
the span of the product states of the form |n,m) with n + m fixed for the individual subspaces. 
Apart from the canonical basis states, the remaining states in the invariant subspaces are entan- 
gled. The states belonging to the invariant subspaces are non-Gaussian states. The know detection 
criteria for entangled states do not identify entanglement in the states belonging to the invari- 
ant subspaces. However, for these non-Gaussian pure states the entanglement criterion is that 

The local modes of sufficient energy are approximate eigenstates of the complete Hamiltonian. 
The entanglement in the evolved state is much larger if the initial state is different from a local 
mode. Though there is a change of quadrature variance during evolution, the states in the invariant 
subspaces do not exhibit single-mode or two-mode quadrature squeezing. The variances are always 
above the minimum allowed limit. The two eigenstates of the total Hamiltonian, belonging to the 
subspace Si correspond to two of the Bell states. Though states of high entanglement cannot be 
generated from the local modes of high energy, the Hamiltonian generates maximally entangled 
states during the evolution of the local modes |0, 1) or |1, 0). 

References 

[1] M. A. Nielsen and I. L. Chuang Quantum Computation and Quantum Information (Cambridge: 
Cambridge University Press) (2000) 

[2] G. Chen et al Quantum Devices Principles, Design, and Analysis (USA: Chapman & 
Hall) (2007) 

[3] C. M. Tesch, L. Kurtz and R. de Vivie-Riedle Chemical Physics Letters 343(2001) 633 
[4] D. Babikov Journal of Chemical Physics 121 (2004) 7577 

[5] A. S. Scott Nonlinear Science: Emergence and Dynamics of Coherent Structures, (London: 
Oxford University Press) (2003) 

[6] A. Chefles and S. M. Barnett Journal of Modern Optics 43 (1996) 709 



19 



L. J. Bernstein Physica D 68 (1993) 174 

S L Braunstein and P van Loock Reviews of Modern Physics 77 (2005) 513 

A. Purusawa et al Science 282 (1998) 706 

E. Ciancio and P. Zanardi Physics Letters A 360 (2006) 49 

J. Perina Quantum Statistics of linear and nonlinear optical phenomena (The Netherlands: 
Kluwer Academic Pubhshers)(1991) 

G. S. Agarwal and R. R. Puri Physical Review A 39 (1989) 2969 

S. M. Jensen IEEE J. Quantum Electron. QE-18 (1982) 1580 

A. Miranowicz and W. Leonski Journal of Physics B: Atomic, Molecular and Optical Physics 
39 (2006) 1683 

W. Leonski and A. Miranowicz Journal of Optics B 6 (2004) S37. 

L. Sanz, R. M. Angelo and K. Furuya Journal of Physics A: Mathematical and General 36 
(2003) 9737 

G. Ariunbold and J. Perina Journal of Modern Optics 48 (2001) 1005 
M. K. Olsen Physical Review A 73 (2006) 053806 

C. Sudheesh, S. Lakshmibala and V. Balakrishnan Journal of Physics B: Atomic, Molecular 
and Optical Physics 39 (2006) 3345 

A. Kowalewska-Kudlaszyk and W. Leonski Physical Review A 73 (2006) 042318 

C. H. Bennett, H. J. Bernstein, S. Popescu and B. Schumacher Physical Review A 53 (1996) 
2046 

M. S. Kim, W. Son, V. Buzek and P. L. Knight Physical Review A 65 (2002) 032323 
Wang Xiang-bin Physical Review A 66 (2002) 024303 

L. Mandel and E. Wolf Optical Coherence and Quantum Optics (Cambridge: Cambridge Uni- 
versity Press) (1998) 

N. Korolkova and J. Perina Optics Communications 136 (1996) 135 

G. C. Gerry and P. L. Knight Introductory Quantum Optics (Cambridge: Cambridge University 
Press) (2004) 



20 



[27] L.-M. Duan et al Physical Review Letters 84 (2000) 2722 

[28] R. Simon Physical Review Letters 84 (2000) 2726 

[29] S. Mancini et al Physical Review Letters 88 (2002) 120401 

[30] E. Shchukin and W. Vogel Physical Review Letters 95 (2005) 230502 

[31] G. S. Agarwal and A. Biswas New Journal of Physics 7 (2005) 211 

[32] H. Nha and J. Kim Physical Review A 74 (2006) 012317 

[33] M. Hillery and M. S. Zubairy Physical Review Letters 96 (2006) 050503 

[34] C. M. Caves Physical Review D 26 (1985) 1817 



21 



